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O ■ ABSTRACT 



The Kennicutt-Schmidt law relates the face-on star formation rate (SFR) per unit 
area with the face-on gaseous column density in nearby galaxies. Applying this relation 
to damped Lya absorption systems (DLAs) of neutral hydrogen column density N > 
1.6 x 10 21 cm" 2 leads to an estimate that three percent of the sky should be covered with 
extended sources brighter than \xy ~ 28.4 mag arcsec" 2 , if DLAs at redshift £=[2.5,3.5] 
undergo in situ star formation. We test this hypothesis by searching the Hubble Ultra 
Deep Field (UDF) F606W image for low surface-brightness features of angular sizes, 
ranging between 0^=0. 25" and 4.0". At z = 2.5 — 3.5, the observed F606W fluxes 
correspond to roughly rest-frame 1500 A and the angular sizes correspond to predicted 
disk diameters of ddia = 2 — 31 kpc. After convolving the F606W image with smoothing 
kernels of angular diameters #kem=#dia> we find the number of detected objects to 
decrease rapidly to zero at #k er n>l "• Our search yields upper limits on the comoving 
SFR densities that are between factors of 30 and 100 lower than predictions, suggesting a 
reduction by more than a factor of 10 in star formation efficiency at z ~ 3. We consider 
several mechanisms that could reduce star formation efficiency at high redshift. We 
find that the cosmological increase with redshift of the critical surface density for the 
Toomre instability may be sufficient to suppress star formation to the levels implied by 



1 Based on observations obtained with the NASA/ESA Hubble Space Telescope. HST is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under NASA contract NASS5-26555. 
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the UDF observations. However, the uncertainties are such that Toomre instabilities 
may still exist. In that case star formation at column densities less than 10 22 cm -2 may 
be suppressed by the low molecular content of the DLA gas. The upper limits on in 
situ star formation reduce the predicted metallicities at z ~ 3 to be significantly lower 
than observed, and reduce the heat input in the gas to be substantially lower than the 
inferred cooling rates. In contrast, the radiative output from compact Lyman Break 
Galaxies (LBGs) with R < 27 is sufficient to balance the comoving cooling rate. This 
leads us to posit that a significant fraction of the DLA population are hosts to more 
compact regions of active star formation, which may be the sources of metal enrichment 
for these DLAs. Such regions are likely to be LBGs. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 

1. INTRODUCTION 

A principal goal of galaxy formation theory is to understand how stars form from gas. This 
condensation process is fundamental for determining the star-formation history of a galaxy, the 
spatial distributions of its stellar populations, and its chemical evolution. Yet it is still not under- 
stood. On small scales, the physical processes are complex and difficult to calculate (e.g. Shu et al. 
1987; Krumholz & McKee 2005), while on cosmological scales, realistic simulations of star forma- 
tion are in addition constrained by limits on numerical resolution (e.g. Cen et al. 2003; Nagamine 
et al. 2004a, 2004b). On the other hand considerable progress has been made in the observational 
sector. Recent multi-color imaging and spectroscopic surveys have succeeded in tracing starburst 
galaxies out to redshifts as large as six (Giavalisco et al. 2004; Bouwens et al. 2004). The majority 
of galaxies found in this way are the Lyman Break Galaxies (i.e. LBGs; e.g. Steidel et al. 2003). 
These are compact (half-light diameters ~ 4 kpc) star-forming galaxies with mean star formation 
rates SFR ~ 40 M Q yr _1 , after extinction corrections are applied (Shaplcy et al. 2003). The LBGs 
contribute a SFR per unit comoving volume p*(z) (Steidel et al. 1999; Giavalisco et al. 2004) that 
in the redshift interval z =[6,2] would consume a mass per unit comoving volume of cold neutral gas 
equivalent to about 10 % of the mass content of visible stars in modern galaxies, i.e., O.lfi^z = 0). 
As a result, reservoirs of neutral gas at z > 2 may be required to fuel these star-forming objects. 

The purpose of this paper is to search for star formation in spatially extended regions that 
serve as such neutral-gas reservoirs. We focus on the damped Lya systems (hereafter DLAs), the 
population of quasar absorption systems with N > 2xl0 20 cm -2 , where N denotes observed H I 
column density (for a review see Wolfe, Gawiser, & Prochaska 2005 [hereafter WGP05]). At z ~ 
3 this column-density threshold guarantees gas neutrality in most cases, which distinguishes DLAs 
from the Lya forest and all other classes of absorption systems in which N < 2xl0 20 cm" 2 and 
the gas is more than 50 % ionized. The neutrality of the gas takes on added significance when it 
is realized that DLAs (1) dominate the neutral-gas content of the Universe at z =[0,5] and (2) at 
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z w 3.5 contain sufficient gas to account for 0.5f2*(z = 0). Furthermore surveys for DLAs reveal 
a large area-covering factor for neutral gas with N > 2xl0 20 cm~ 2 (Wolfe et al. 1995; Prochaska, 
Herbert-Fort, & Wolfe 2005 [PHW05]). In the redshift interval z= [2.5,3.5], a redshift range in 
which the SFR history of galaxies is well determined, DLAs cover one third of the sky, which leads 
to a striking conclusion: the sky should be "lit up" with star formation. Specifically, a subset of 
DLAs with N > 1.6xl0 21 cm -2 should cover three percent of the sky and have surface brightnesses 
brighter than 28.4 mag arcsec -2 (see below). This follows if one presumes the Kennicutt-Schmidt 
law (Schmidt 1959; Kennicutt 1998a, b) holds at high redshifts. In that case star formation occurs in 
the presence of cold atomic and/or molecular gas with a SFR per unit area projected perpendicular 
to the disk given by 



where N± is the H I column density perpendicular to the disk. In nearby galaxies, K=KKenn 
=(2.5±0.5)xl0" 4 M yr" 1 kpc" 2 , /3=1.4±0.15, and the scale factor JV c =1.25x 10 20 cm~ 2 (Ken- 
nicutt 1998a,b). The threshold column density Nf lt is observed to range between 5xl0 20 cm~ 2 
and 2xl0 21 cm~ 2 (Kennicutt 1998b) and is usually associated with the threshold condition for the 
Toomre instability. We ignore the molecular component of the gas owing to the very low H2 area 
covering factor of DLAs (e.g. Ledoux et al. 2003). 

In this paper, we test the hypothesis that star formation proceeds throughout the absorbing 
gas at the projected rates given by Eq. 1. That is, we consider whether star formation at high z 
occurs in sites other than compact LBGs. In § 2 we apply Eq. 1 to high-z DLAs and compute their 
expected surface brightnesses. We compute for the first time the SFR per unit comoving volume 
for DLAs modeled as randomly oriented disks. We then calculate the number of DLAs expected 
to occupy the Hubble Ultra Deep Field (UDF; S. V. W. Beckwith et al. 2004, in preparation). In 
§ 3 we describe a search in the UDF for extended regions with low surface brightnesses. We use a 
matched kernel technique optimized for detecting faint emission from extended objects and discuss 
the results of our search in § 4. In § 5 we describe the implications of these results. A summary 
and concluding remarks are given in § 6. 

Throughout this paper we adopt a cosmology with (Om, ^a, h)= (0.3, 0.7. 0.7) (Spergel et al. 
2003). We adopt the AB magnitude system and refer to magnitudes and surface brightnesses 
deduced from the F606W image with the ACS camera on HST as follows: V = AB(F606W) and 

fly = MF606W- 




(1) 



2. 



THEORETICAL FRAMEWORK 



Before describing the observational results and analysis we discuss observational properties 
predicted for DLAs relevant for detecting them in emission. 
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2.1. Comoving SFR Densities Predicted by the Kennicutt Schmidt Law 

In order to estimate the surface brightnesses of DLAs with projected SFRs predicted by the 
Kennicutt-Schmidt law, we assume they are disks inclined to the plane of the sky by inclination 
angles i. In that case the observed intensity at frequency uq is given by 



47r(l + z) 3 cos(i) 



^(o = ,:tt„,* > ( 2 ) 



where (£^)j_ is the luminosity per unit frequency interval per unit area projected perpendicular 
to the plane of the disk and v = (1 + z)vq\ owing to the low dust-to-gas ratios of most DLAs 
(Pettini 2004), we ignore the effects of extinction. In the case of FUV radiation (with A « 1500 A), 
(£„)±=C where C=8xl0 27 ergs s" 1 Hz-^M© yr- 1 )- 1 /(3.08xl0 21 cm/kpc) 2 = 8.4xl0" 16 

ergs cm -2 s _1 Hz _1 (Mq yr" 1 kpc -2 ) -1 , where the calibration is insensitive to wavelength in the 
FUV portion of the spectrum (Madau, Pozzetti, & Dickinson 1998). Because N±=cos(i)N, we find 
that, for a fixed observed column density, I Vo (i) averaged over 0° < i < 90° is given by 

< ^ >= M if z) 3p A*=K{N/N c f . (3) 

Note that tp* is not the average SFR per unit area projected along the line of sight, <(^*)j_/cos(z)>. 
Rather it is an effective projected SFR relating the observed column density, N, to the surface 
brightness of the DLA. As an example, consider a DLA with iV =1.6xl0 21 cm" 2 , which is about 
twice the mean value of N for the statistical sample of over 625 DLAs (PHW05). From Eq. 3 
we find V ; *=8.9xl0 -3 M yr" 1 kpc" 2 (about twice the local rate in the Galaxy). At z = 3 this 
corresponds to an AB surface brightness in the V band, nv= 28.4 mag arcsec -2 . Here we assumed 
the values for K, /?, and N c cited above (Kennicutt 1998a,b). While objects this faint are beyond the 
sensitivity of the Hubble Deep Field (HDF) and the GOODs survey, they are within the sensitivity 
of images acquired by the UDF (Bouwens et al. 2004). 

We next compute the comoving SFR density p*(z) predicted by the Kennicutt-Schmidt law in 
order to compare it with the empirical values determined from the UDF (see § 3). The calculation 
was originally carried out by Lanzetta et al. (2002: see also Hopkins et al. 2005). Here for the first 
time we account for inclination effects that must be present if, as we assume, DLAs are disk-like 
structures or any type of gaseous configurations with preferred planes of symmetry, such as those 
predicted in high-resolution numerical simulations (e.g. Razoumov et al. 2005; Kravtsov 2003). 
For disks with inclination angles between % and i + di the comoving SFR density due to the intrinsic 
differential area dA± of the disk is given by 

d 2 p* = n co dA±cos(i)x [(ifj*)±/cos(i)]xsm(i)di , (4) 
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where n co is the comoving density of disks. Noting that the intrinsic column-density distribution 
of the disk is defined by the relation 

g(N ± ,X)dN ± = (c/H ) nco (X)dA ± , (5) 
and cos(i)=N ± /N and sm(i)di=(N ± /N 2 )dN we find that 

rN max rmin(N ,N) 

p4>N,X) = (H /c) dNK(N/N c f dN ± g(N ± ,X)(Nl/N 3 )(N ± /N)^ 1 , (6) 

Jn JN min 

where g(N±,X) and the observed column-density distribution function f(N,X) are related by 

/min(A r o,A r ) 
dN ± g(N ± ,X)(Nl/N 3 ) (7) 

(Fall & Pei 1993; Wolfe et al. 1995) and X(z) is the absorption distance (Bahcall & Peebles 1969). 
In deriving these equations we implicitly assumed N±_ to be a monotonically decreasing function 
of radius with a maximum value of Nq. Note, in the spherically symmetric limit the (N^/N) 13 ' 1 
term in Eq. 6 would be replaced by 1 and consequently the expression for p*(>N, X) would be the 
same as the Lanzetta et al. (2002) result (see also Hopkins et al. 2005). 

Equations 6 and 7 show that one must obtain g(N±,X) from f(N, X) to compute p\(>N, X). 
The double power-law fit to the SDSS data resulted in f(N, X) =k 3 (N/N d ) a where k 3 =(1.48±0.07) x 10" 
cm 2 , a=a 3 =-2.00±0.06 for N < N d and a=a A =-6.00^;^ at N > N d , where JV d =(3.16±g;||) x 10 21 
cm" 2 (PHW05). It is natural to equate N d with Nq in which case Eq. 7 shows that projection ef- 
fects alone result in a^= —3.0. This solution is clearly consistent with the data, given the large 
uncertainties in «4 (note, p\[> N, X] is independent of a^). As a result this fit leads to g(N±, X)= 
k3(N_i_/N )~ 2 at N± < N and g(N±,X)=0 for N± > N . One could envisage an alternative 
solution in which A*o ~ 10 22 cm" 2 and in which dlogg(Nj_, X)/dlogN± decreases below — 3.0 as 
A^j_ increases toward A^o- This is observed locally and is presumably due to the conversion of H I 
into H2 at large column densities (Zwaan et al. 2006). However, the molecular fraction does not 
increase with A^ in DLAs (Ledoux et al. 2003), which is likely due to low dust content. Since the 
gamma distribution fit, f(N,X) oc (N/N^expi-N/N^) (PHW05), leads to negative g(N ± ,X) 
for N± > (3+a2)A 7 , we conclude that our double power-law fit to f(N,X), with a 4 = — 3.0, is 
the most plausible solution for DLAs, provided they are disk- like structures, and we adopt it here 
and in what follows. 

The resulting cumulative comoving SFR density is depicted by the blue curve in Fig. 1. We 
plot p*(>ip*,X) for the range in ip* corresponding to N = [N m i n , N max ], where here and unless 
otherwise noted, A r m j n =2xl0 20 cm" 2 and N max =10 22 cm -2 . By adopting a value for N m i n that is 
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Fig. 1. — Cumulative comoving SFR density versus "0*. Curve is /0*(>V'*) predicted for Kennicutt- 
Schmidt law with K=KKenn- The range of ip* corresponds to 7V=[2xl0 20 ,lxl0 22 ]. 
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lower than the range of threshold column densities Nf lt observed for nearby galaxies (Kennicutt 
1998b), we overestimate p*(>N m i n ). However, the upper limits on p\ estimated in § 3.3 are valid 
only for N higher than the local values of N^ lt , and as a result comparison between theory and 
observation occurs at column densities where the Kennicutt- Schmidt law is well established. We 
also assume z=3 and use Eq. 3 and the Kennicutt parameters to convert N to tp*. In agreement 
with Hopkins et al. (2005) we find that p* inferred for DLAs is predicted to be somewhat lower 
than deduced for LBGs (Steidel et al. 1999; Giavalisco et al. 2004), or for DLAs from the C II* 
technique (Wolfe et al. 2003a [hereafter WGP03]). 



2.2. Expected Number of DLAs in the UDF 

How many DLAs are expected to populate the UDF? The answer depends on three factors: (1) 
the column-density distribution function f(N,X), (2) the redshift search interval z = [zi ow , z^igh], 
and (3) the linear sizes of the DLAs. Suppose all DLAs were equal in size. Because angular 
diameter is an insensitive function of z at z > 1.2, we assume that all DLAs in the search volume 
have identical angular diameters, Odla- In that case the number of DLAs with column-densities 
greater than N is given by 

-^j-A / dX dN'f(N',X) , (8) 

™dla/ a J Jx(z low ) Jn 

where ujjjdf is the solid angle subtended by the UDF. The results shown in Fig. 2 were obtained as- 
suming loudf =10 arcmin 2 , zi ow =2.5, and 2^^=3.5. The four curves were computed for 9dla=0-5" , 
1" , 2", and 4", which correspond to linear diameters of 4, 8, 15, and 31 kpc at z=3 that span the 
range of sizes predicted for DLAs, which can reach diameters as low as 2 kpc (e.g. Haehnelt et al. 
1998; Nagamine et al. 2005; Razoumov et al. 2005; Prochaska & Wolfe 1997). To evaluate the 
integral in Eq. 8 we combined the redshift evolution of dAT/dX with the A-dependence of f(N, X) 
determined from the full sample of DLAs found in the SDSS survey (PHW05). The figure shows 
the UDF should be populated by large numbers of DLAs. In our example of N >1.6xl0 21 cm -2 
we find that over 100 DLAs would be present with py < 28.4 mag arcsec -2 if Odla = 4". Of 
course, not all DLAs will have the same size, and most should have Odla < 4", which results in 
an increase in M. The point is that the UDF should contain hundreds to tens of thousands of 
high column-density DLAs which according to Eq. 3 should be detectable. By comparison, were 
the LBG luminosity function (Adelberger & Steidel 2000) extrapolated to an approximate UDF 
point-source limit of U=30.5, the UDF would be populated with 1800 LBGs in the redshift interval 
z=[2.5,3.5]. 
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Fig. 2. — Plot showing number of DLAs with column densities exceeding N predicted to populate 
the UDF in the redshift interval z= [2.5,3.5]. Magenta, blue, red, and green curves are for DLAs 
with 9 DLA =0.5", 1", 2", and 4 ". 
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3. ANALYSIS 

We now describe the results of our search for emission from DLAs. 

3.1. Object Detection 

To determine p*(>ip*) (here and in what follows we let p*(>?p*)=p*[>ip*,X(z = 3)]) em- 
pirically, we searched for extended, low surface-brightness emission in the UDF. The search was 
performed with images acquired with the F606W filter. The wavelength centroid of the F606W fil- 
ter approximately matches the rest-frame FUV wavelength of 1500 A for the redshift search interval 
z= [2.5,3.5]. This is the standard rest-frame wavelength for which the LBG luminosity function and 
comoving SFR densities have been determined (Steidel et al. 1999; Adelberger & Steidel 2000). 
Furthermore, this redshift interval contains the largest sample of known LBGs, which we shall use 
for comparison with our results, and the F606W image is the most sensitive one obtained in the 
UDF. Because the UDF does not have the i7-band sensitivity required for this project and the 
Lyman limit discontinuities of z ~ 3 objects occur further blueward of the blue filter response, pho- 
tometric redshifts for objects in this redshift range are subject to catastrophic errors. We therefore 
apply the data only for determining upper limits on p*(>V>*) by assuming that all objects lie in this 
redshift range. 

We performed object detection in the drizzle-combined F606W image using SExtractor (Bertin 
& Arnouts 1996). The object detection criteria are set according to the following steps. First, we 
require the presence of flux over a minimum, contiguous area comparable to the size of an adopted 
kernel. Next, we adjusted the detection threshold to the lowest value where zero detections are 
found in the negative image. This procedure was first performed in the original F606W image with 
a smoothing kernel matching the PSF of the image. Specifically, we let the FWHM of the kernel 
#kern=#PSF = 0.09", in order to identify all the high surface brightness (HSB) objects; i.e., objects 
with ny < 26 mag arcsec -2 . We then masked out all these known objects in the image by setting 
the associated pixels to the median sky value, in order to produce a masked-F606W image. To 
search for extended low surface brightness objects, we smoothed the masked-F606W image using 
circular Gaussian kernels with #kcrn=#dia where #dia is the FWHM of the DLA image and let the 
aperture diameter equal #kcrn- The advantage of smoothing is that in the limit of random noise 
the SNR of the surface brightness of the smoothed image is enhanced relative to the unsmoothed 
image by the factor #kern/#PSF (e.g. Gonzalez et al. 2001; Irwin et al. 1985). We emphasize that 
removing HSB objects detected in the unsmoothed image from all subsequent searches allows us 
to search for extended low surface-brightness features closer to these HSB objects in smoothed 
images. Consequently we are testing the hypothesis that emission from DLAs arises only from star 
formation in the gas detected in absorption, since the Kennicutt-Schmidt law predicts py to be 
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fainter than 26 mag arcsec 2 for N=[N m - m , N max ] and z = [2.5,3.5]. 2 

At z >2 most HSB objects are LBGs. While the brighter, V < 25, LBGs have half-light radii, 
r/j/«2kpc (Giavalisco et al. 1996), recent studies with the UDF indicate is smaller for fainter 
objects and that rhi decreases with increasing redshift. Specifically Bouwens et al. (2004) find that 
the mean half-light radius, < rhi >=0.9 kpc for their sample of intermediate magnitude, Z850,AB < 
27.5, objects. We shall assume a conservative upper limit < r^i ><1.5 kpc for LBGs in the UDF 
with V < 30. The light profiles of LBGs are well fitted with exponentials. Since exponentials have 
FWHM linear diameters equal to 2rw/2.5, the expected FWHM angular diameters for most LBGs in 
our sample, 9lbg<®-13" . By contrast the angular diameters predicted for the light-emitting regions 
of DLAs are systematically larger. Although the observational evidence for redshifts exceeding 1.9 
suggests 2.0 kpc < ddi a < 15 kpc (where <idia is the FWHM linear diameter), the data sample 
contains only three objects (WGP05). On the other hand, theoretical arguments suggest DLAs 
are larger than LBGs. While the predicted sizes of DLAs are smallest for CDM models, most of 
the predictions for z ~ 3 indicate d^u > 1-9 kpc. Comparison with the numerical simulations of 
Nagamine et al. (2006) shows that 70 % is a conservative lower limit for the fraction of z=3 DLAs 
with c?dia > 1-9 kpc. Similarly, Haehnelt et al. (2000) predict at least 80 % of these DLAs to 
have ddia > 1-9 kpc, while the lower limit computed by Mo et al. (1998) is 95 %. Because other 
models (e.g. Prochaska & Wolfe 1997; Boisser etal 2003) predict even larger sizes, current ideas 
about DLAs suggest that the bulk of DLAs predicted in current models have linear diameters larger 
than 1.9 kpc, which corresponds to 0.25 " for the redshift interval z=[2.5.3.5]. We shall revisit the 
possible connection between DLAs and HSB objects in § 5.2. 



3.2. Search Results 

Using the #kern=0.09 " smoothing kernel, we identified roughly 11,000 objects brighter than the 
survey limit of V < 30.5 and with (nv) < 26.6 mag arcsec -2 in the original F606W image over the 
central 10 arcmin 2 of the UDF. Interestingly, none of the objects satisfied our search criteria for in 
situ star formation throughout the DLA gas; i.e., low surface brightness, [Xy > 26 mag arcsec -2 , and 
angular diameter exceeding 0.25 ", which corresponds to the smallest linear diameter 3 predicted 
for a substantial fraction of model DLAs (§3.1). We considered the "tadpole galaxies" identified at 
z < 4 in the UDF by Straughn et al. (2006) as potential candidates because they typically consist 
of multiple compact components distributed over an area ~ 0.1x1 arcsec 2 : we confirm the presence 
and properties of all 163 "tadpole galaxies" identified by Straughn et al. (2006) and measure their 



2 We chose this surface-brightness limit because at z = 3 it corresponds to N < 8xl0 21 cm 2 , the highest column 
density yet measured in the SDSS DLA survey (PHW05). 

3 Most discussions about DLA sizes refer to impact parameter 6 rather than diameter d. In the simple case of 
randomly oriented circular disks, the median impact parameter b=d/2\/2, which corresponds to b — 0.67 kpc for #kem 
=0.25 ". 
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average V surface brightness < fly > = 26.6 mag arcsec~ 2 . These authors detected between 50 
and 60 galaxies with provisional redshifts between 2.5 and 3.5, which leads to an area covering 
factor /a ~ 10~ 4 . By contrast /^=0.33 for all DLAs in this redshift interval, and //i=0.13 for 
DLAs with TV > 5xl0 20 cm -2 , the H I column density above which the Kennicutt-Schmidt law is 
firmly established (Kennicutt 1998a,b). As a result the "tadpole galaxies" cannot be the sites of in 
situ star formation in most DLAs. On the other hand, the relatively high surface brightness of the 
compact components in these galaxies is consistent with star formation occurring in regions with 
N > 5xl0 21 cm -2 that follow the Kennicutt-Schmidt law. 

While we have not attempted to measure the redshifts of the remaining high surface-brightness 
sources in our sample, the majority of objects with fiy < 26.0 mag arcsec -2 are likely to be faint 
blue galaxies with z < 1.5. This conclusion is based on the luminosity function presented by 
Blanton et al. (2005), which includes a double Schechter function to account for the steep rise in 
galaxy density at the faint end. At r' = 28 about 10 % will be at z > 3, with a median redshift of 
z = 1, while at r'=30 about 30% will be at z >3 with a median redshift of 1.5. These numbers are 
relevant to the F606W image because its wavelength centroid is similar to that of the r' filter. 

To increase the probability for detecting extended low surface brightness objects, we searched 
the masked F606W image by convolving it with smoothing kernels matched to a range of predicted 
0dia's, from ^dia=0.25" to #dia=4.0" (see § 2.2). The number of objects found in the smoothed 
images versus the smoothing kernel size is presented in Fig. 3. The figure shows that the number 
of detected objects decreases rapidly with increasing 6>ke m : while 166 objects are detected when 
#kern=0.25", only one object is found when 0kern=LO", and no objects are found for 9^ CTn > 1.(1' . 
In § 3.3 we shall argue that the steep fall off with #kem indicates the presence of compact sources 
with light profiles that decline steeply with radius. 

The implication is that spatially extended sources of low surface-brightness emission are rare in 
the UDF. Fig. 4 displays multi-color images of the 12 objects detected for 0k ern =O.5". Comparison 
between the smoothed and unsmoothed F606W image demonstrates the efficiency of the matched 
kernel technique for detecting extended objects that are exceedingly faint. The detection of objects 
with 6*kern = 0.5" supports the conclusion that had they been present, faint objects with #dia > 1" 
would also have been detected with kernels appropriately matched to #dia- Lacking photometric 
redshifts, one cannot determine the nature of the 12 objects but it is likely that they are low 
surface-brightness dwarf galaxies at redshifts z ~ 1. The blue colors, central surface brightnesses, 
and linear diameters of the sample compiled by van Zee et al. (1997) would give rise to V band 
magnitudes, surface brightnesses, and angular diameters similar to those detected here, if the 
sample were redshifted to z « 1. Because of the similarity in V magnitudes, we also considered 
whether our detected objects belong to the class of IRAC selected galaxies found in the UDF by 
Yan et al. (2002). But this is unlikely since the IRAC galaxies are brighter at Zs5o,ab magnitude 
than at V, whereas all of the objects that we detected in V with #kern=0.25 " to 1.0 " are either 
fainter or undetected in Z850,AB- As a result, the density obtained by assuming these objects are 
in the redshift interval z= [2.5,3.5] acts as a conservative upper limit to the density of DLAs with 
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Fig. 3. — Histogram showing number of recovered DLAs as a function of #kem- Steep decrease 
with increasing #kern indicates that increase in aperture size results in increase in sky noise but 
not in signal. This implies the predominance of compact galaxy light profiles reminiscent of dwarf 
galaxies at intermediate redshift. Part of fall-off is also due to the observed intrinsic decrease in 
galaxy counts with increasing flux. 



-13- 



^dia=0-5" in the same redshift interval. Similar conclusions hold for all DLAs with 6*di a < !"• 

Therefore, in the cases of null detections, #dia > 1 "> we use cumulative Poisson probabilities 
to place a 95% confidence upper limit on N of A/95 = 3 for the number of low surface-brightness 
DLAs populating the UDF in the redshift interval z = [2.5,3.5]. In the case of a single detection, 
#dia = 1 "j we place a corresponding 95 % confidence upper limit of A/95 = 4.7. For the multiple 
detections at #dia < 1 " ■> the values of A/95 are given in Table 1. We emphasize that the upper limits 
on N are valid only for objects brighter than some threshold magnitude and surface brightness, 
which we determined by a Monte-Carlo technique described in the following subsection (§ 3.3). 

3.3. Search Completeness 

To determine the sensitivity of our searches for extended low surface-brightness objects, we 
performed a series of Monte-Carlo simulations. First, we generated test objects of different bright- 
nesses for exponential light profiles of different scale-lengths. The angular exponential scale-lengths 
r s were adjusted such that the FWHM (=[21n 2]r s ) = #dia- Next, we placed the test objects of 
a given magnitude and scale-length at 1000 random positions in the masked-F606W image, and 
smoothed the images with kernels in which #k ern = #dia in the optimal case of matched kernels. 
Finally, we performed the detection procedure following the same detection algorithm described 
in § 3.1 for finding low surface brightness objects. We also repeated the simulations in which the 
synthetic galaxy light profiles were elliptical exponentials with 6/a=0.5 and with position angles 
varying between 0° and 90°: the difference with results for circular apertures was insignificant. 

We measured the photometric properties of the recovered galaxies and computed the recovery 
fraction, e. This procedure was repeated by increasing the input magnitude V m until the recovered 
number of galaxies decreased to some critical value. The results of this procedure are illustrated in 
Fig. 5, which plots e versus V m . In principle one could determine the threshold value of V by raising 
V m until only A/95 objects are left in the image. For those cases where A/95 << 200 (i.e. e << 0.2) 
we found that the measured magnitudes, V meas , were sufficiently faint that overlap between the 
synthetic galaxies and the low-surface brightness wings extending from the masked real galaxies 
led to spurious over-estimates in flux determinations for the synthetic galaxies (in most cases this 
effect results in y meas to be brighter than V m ). This is evident in the saturation of the e versus 
V m curves in Fig. 5: the flattening at e < 0.2 for #dia=4 " and at e < 0.1 for the other values of 
#dia illustrates this effect. On the other hand when e > 0.2, the V m were sufficiently bright that 
the output V meas were reliable. As a result we assume the threshold magnitude is given by V m for 
e = 0.2 and the threshold surface brightness is given by the corresponding value of \i l y (hereafter 

Fig. 4a. — multi-band images of galaxies detected with 0.5 " kernel. Column 1 shows smoothed 
F606 W image, col. 2 shows F606W masked image. Remaining columns show unsmoothed images 
for designated filters. Images are 5 " on each side. 
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Fig. 4b. — Same as 4a 



Fig. 4c. — Same as 4a 
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Table 1: Results of UDF Search in F606W Image 



"Light profile exponential scale-length for synthetic galaxy used to compute threshold value of ip* for all cases and 
threshold p* in the case of null detections. 
6 FWHM of Gaussian smoothing kernel 

c Number of objects detected in the UDF with given value of 0kem 

d 95 % confidence upper limit on number of objects in UDF with given #kern 

e Threshold input V magnitude, when e=0.2. Used to compute upper limits on SFRs for kernels with null detections. 
Not used for kernels with positive detections. 
* Threshold input V band surface brightness when e=0.2 

9 Upper limit on SFR. For A/det = this corresponds to threshold V magnitude, and for Adet =fc this corresponds the 

SFR averaged over the A/det galaxies. 
''Lower limit on tp* (z=3.5) corresponding to threshold fj,v 
1 95 % confidence upper limit on p* 
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Fig. 5. — Recovery fraction e versus input magnitude V m for increasing values of 9^ ern for the 
matched kernel condition, #kcrn=$dia- Notice how V' m at a given recovery fraction, e=0.2, decreases 
with increasing #kern- This results from noise increase due to increase in aperture size. The flattening 
observed at e=0.2 is due overlapping of synthetic galaxy profiles with low brightness wings extending 
from masked galaxies. 
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we refer to these as V thresh and p l y resh ). We emphasize that the e=0.2 criterion is conservative 
as it requires the detection of 200 objects, which exceeds A/95 f° r a U values of O^u- The resulting 
threshold magnitudes and surface brightnesses corresponding to the values of #kem used for galaxy 
detection are shown in Table 1. 

A possible concern over this procedure arises from our use of 10 3 synthetic galaxies even though 
Fig. 2 shows that A/", the number of DLAs predicted for the UDF, can be as large as 10 5 . However, 
the threshold values of V m and p,™ inferred from the simulations are independent of M . The reason 
is that confusion noise due to fluctuations in the number of synthetic galaxies per aperture beam 
does not increase with increasing N: for DLAs with a fixed area covering factor /^(DLA) the 
number of galaxies per "beam" area, vr6'] 5 - CI . n 2 /4, equals fA{DLA){9^ eTn /dg[ a ) 2 , which is independent 
of N and is less than 1 for the case of matched kernels assumed here. Other systematic errors such 
as that due to overlap between the galaxy light profiles is negligible because the area covering factor 
of 200 synthetic galaxies is less than 6 %. Therefore, the principal systematic error in determining 
the thresholds is caused by overlap between the synthetic light profiles and the faint wings of 
masked real galaxies (see above). 

We next used the simulations to interpret the steep fall-off in the number of detected galaxies 
with 6>kern (Fig. 3). Part of the fall-off is due to the loss of sensitivity resulting from the increase 
of #kcm above #dia- This is illustrated in Fig. 6, which shows that V thresh decreases with increasing 
$kern f° r $kcm > #dia- Therefore, a loss of sources will result from a decrease in survey depth with 
increasing #kem- A fall off of sources with decreasing y thresh can a l so be attributed to the intrinsic 
shape of the source counts-magnitude relation. From known field number counts (Metcalfe et al. 
20001), we estimate that for every half magnitude deeper, we gain about 20% more objects at V 
~28. If we adopt this rate and our survey depth estimated from an exponential profile, which went 
from V thresh =29.0 to V thresh =28.7 when we increase the smoothing kernel size from 0.3" to 0.4", 
we estimate a loss of only about 13 % of the targets if the decline is due to survey depth alone. 
Because Fig. 3 indicates a more significant drop than this, we suggest within one exponential scale 
length the intrinsic light profiles of the detected galaxies are steeper than an exponential function, 
which is consistent with the dwarf galaxy hypothesis. 



3.4. Estimated Comoving SFR Densities 

We computed empirical upper limits on p*(>ip*) in the following manner. First, in the case of 
null detections, #kern=2" and 4" , we used the threshold magnitudes y thresh in Table 1 to determine 
the minimum luminosity, J J r ^ %n detectable out to z^gh — 3. 5. We then converted L™ n into a SFR 
by assuming SFR=1.25x lQ- 28 L™ in (Madau et al. 1998), and let p\=J\f 95 xSFR/ AV co to find the 
95% confidence upper limit on p\, where AV^, the comoving volume of the UDF, equals 3.3 xlO 4 
Mpc 3 for the redshift interval z= [2.5,3.5]. We then used Eq. 3 to compute the threshold value of 
V>* from the threshold surface brightnesses py resh in Table 1. The resulting p*(>ip*) are plotted in 
Fig. 7 and entered in Table 2 for both values of #kcm- In the case of positive detections, 0kem < 1", 
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Fig. 6.— Threshold magnitude V thresh versus 6» kcrn for fixed DLA diameters, d i a = 0.25", 0.5", 
and 1.0 ". Decrease of V thresh with increasing 6>kcrn at #kcm > #dia due to increase in sky noise in 
larger apertures. Decrease of V thresh with decrease in #kem at ^kem < #dia due to effect of smoothing 
kernel alone. 
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Fig. 7. — Cumulative comoving SFR densities versus Green data points correspond to upper 
limits derived for 6> kcrn =4", 2", 1", 0.6", 0.5", 0.4", 0.3", and 0.25 " in counter-clockwise order 
starting with (logi/>*,log/j*)=(— 2.41, —2.97): matched kernels assumed in deriving SFR thresholds. 
Blue (starred) data points correspond to unmatched kernel with 6^=1.0 " and #kcm = 

0.5" and 

2.0" in order of increasing ip*. Blue curve is p*(>ip*) predicted for Kennicutt- Schmidt law with 
K=Kxenn- Solid Red curves are for K=0.1KKenn and K=0.0lKKenn- The range of ip* for all 
solid curves corresponds to iV=[2xl0 20 ,lxl0 22 ]. Dot-dashed red curve is same as solid red curve 
for K=0.1KKenn except that N max =2x 10 22 cm -2 . 
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we assumed all the detected objects were at z=3 (since the average redshift of all of the detected 
objects is more likely to equal the average redshift rather than the maximum redshift, z^igh, of 
the UDF search volume) and added the inferred values of L v to obtain the total SFR. Because 
A/95 exceeds the number of detected galaxies, A/det> the luminosities of the A/95 — Adet galaxies are 
unknown. We addressed this issue by assigning the mean SFR of the detected sample <SFR> to 
each excess galaxy. Therefore, in these cases p*= A/g5<SFR>/AV r co . Note, whereas the threshold 
values of ^* were computed with the same method used in the case of null detections, in this case 
the limits on p\ depend on the V magnitudes measured for each detected galaxy rather than on 
the value of V thresh used in the case of null detections. 

We wish to comment on two more points about Fig. 7. First, the horizontal green arrows 
directed toward increasing tp* emphasize that each upper limit on p*(>Y>*) is valid for projected 
SFRs > Second, the upper limits on p\ decrease with decreasing #kcrn between #kcm = 4.0 " 
and 1.0 " because A/95 is essentially unchanged while the sky noise contribution to the measured 
flux is diminished in the smaller kernels. This is reversed to an increase in p\ with decreasing #kem 
for #kem < 1-0 " because the rapid rise in A/95 with decreasing #kem (Fig.3) offsets the reduction 
in sky noise. The corresponding ip* coordinate increases with decreasing #kem i n every case due to 
the increased error in surface brightness with decreasing #kcm- 



4. DISCUSSION OF RESULTS 

The most striking feature of Fig. 7 is that the comoving SFR densities predicted by the local 
Kennicutt-Schmidt law (the blue curve) exceed the 95 % confidence upper limits established for 
the UDF by at least a factor of 30 for #kem = 0.25 ", 0.30 ", and 4.0 ", and by more than a factor 
of 100 for other values of #kern- 

We believe these discrepancies are real for the following reasons. 

(1) the discrepancy is established for a wide range of #kcrn values, indicating the results are 
not overly sensitive to the values of #dia- 

(2) The upper limits are likely to be conservative; i.e., the locations of the data points are 
likely to occur at values of V>* and p\ lower than in Fig. 7. In the case of null detections, 6>kern = 2.0 
" and 4.0 the location of the upper limits in the ip*, p\ plane were determined by the threshold 
surface brightnesses and magnitudes determined in our simulations. In § 3.3 we discussed why our 
recovery criterion of 200 objects led to large values of ip* and p\. In the case of positive detections, 
#kern < 1-0 the threshold values of ip* are again determined by the simulations, but the upper 
limits on p\ are set by the magnitude determinations of the detected objects. Because a significant 
fraction of these objects are likely to be intermediate-redshift dwarf galaxies, the contribution of 
all the detections to the total SFR in AV co is a conservative upper limit. 

(3) The discrepancy is not overly sensitive to the matched-kernel assumption described in § 3. 
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The two blue stars in Fig. 3 correspond to tp* and p* inferred for #dia=l " and #kem = 0-5 " and 
2.0 " . The small changes in threshold ^* support our conclusion that the discrepancy between 
empirical and predicted estimates in the ip*, p* plane is robust. 

(4) Each of the upper limits on p* in (Fig. 7) applies to DLAs of the same size. But in 
a realistic scenario, DLAs contributing to p* would span some range in size (or in mass since 
size is correlated with dark- matter mass in most models). Springel & Hernquist (2003) found 
that p*=^MPcritf^° S(M, z)dM where p cr u is the current critical cosmological density and the 
multiplicity function S(M, z) is the product of the comoving density of dark-matter halos and the 
SFR per unit halo mass. Therefore, the SFR density predicted for a given smoothing kernel is 
given by (p;) 9kcrn = [J™+ S(M, z)dM / J °° S(M, z)dM]p\ where M_ and M + are the mass limits 
corresponding to the range in #dia detectable with a given #kcm- As a result the only relevant upper 
limits are those corresponding to kernels sensitive to masses contributing the bulk of the predicted 
p'*; i.e., kernels for which (p*)e kcrn is not much less than p*. Our simulations indicate that the 
smoothing kernels maintain maximal sensitivity for DLAs with #dia=(0.5 — > 1.5)0kcm- Thus, the 
kernel with 0kern=2" is sensitive to z=3 DLAs with c^ a =8kpc to 24 kpc. In the spherical collapse 
model the corresponding range in dark-matter masses is 1O 11 M0 to 1O 13 M0 if we assume the virial 
radius r2oo ~ 10(d^ o /2) (see § 5.1). Fig. 9 in Springel & Hernquist indicates that halos in this mass 
range contribute about 50 % of the predicted p* at z=3. This model would be ruled out by our 
data since the predicted p* is 100 times larger than the empirical upper limit for 0kem=2". While 
more rigorous methods are needed to verify this result, our preliminary analysis indicates that the 
results presented here are valid in the case of distributed DLA sizes. 

(5) The values of p'*(>Y>*) predicted for the local Kennicutt parameters, K and (3, are likely to 
be higher than depicted by the blue curve in Fig. 7. While the predicted p'*(>Y>*) is well determined 
by the accurately measured parameters such as the slope 03, normalization k^, and turn-over column 
density Nj, it also depends on the uncertain value of N max . We chose N max =lxW 22 cm -2 because 
the largest measured column density in the SDSS sample is A=8xl0 21 cm~ 2 (PHW05). The blue 
curve is a lower limit because p*(>ip*) would increase if N max were raised above lxlO 22 cm~ 2 , 
which cannot be ruled out. Furthermore the predicted p*(>ip*) would also increase if DLAs arc 
not the planar objects assumed in § 2, 

(6) A potential problem with our results is that while the UDF measurements sample scales 
exceeding 1 kpc, the expression for p\ in Eq. 6 depends on f(N,X) which is based on absorption- 
line measurements that sample scales ~ 1 pc (Lanzetta et al. 2002). The difficulty is that the 
Kennicutt-Schmidt law has been established only on scales exceeding 0.3 kpc (Kennicutt et al. 
2005). However, because f(N,X) typically depends on over 50 measurements per column-density 
bin of width AlogA=0.1 (PHW05), Eq. 5 implies that the corresponding comoving SFR density Ap* 
is a statistical average over the differential area AA^(Ho/c)f(N,X)AN/n co (assuming spherical 
symmetry). We find that AA m 0.5(Ad/A)(# d ia/l arcsec) 2 kpc 2 . As a result AA w 2 kpc 2 for 
#dia=l " at the median DLA column density A=8xl0 20 cm~ 2 . Since our calculation corresponds 
to the cumulative quantity p*(>N), the effective statistical area probed by our expression will 
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exceed a few kpc 2 . Therefore, while individual measurements sample scales less than 1 pc, p* is 
statistically established on scales for which the Kennicutt-Schmidt law has been established at z=0. 

5. IMPLICATIONS 

The discrepancies with predictions by the Kennicutt-Schmidt law have implications for the 
efficiency of star formation in DLAs. They also have broader implications concerning metal pro- 
duction and energy balance in DLAs and their relationship to LBGs. We discuss each of these 
topics in turn. 

5.1. Star Formation Efficiency in DLAs 

The discrepancy between predicted and measured values of p\ suggests that star formation is 
less efficient in DLAs; i.e., for a given N, ijj* in DLAs is lower than in nearby galaxies. This point 
is illustrated by the solid red curves in Fig. 7, for which the Kennicutt parameter K is reduced 
by factors of 10 and 100 below the local value, K=K-^ cmi . As K decreases, the solutions shift 
toward decreasing ip* and p* along a 45° diagonal in the (log^*, log/3*) plane as both ip* and p\ 
are proportional to K. Because the upper limits on p\ apply only to projected SFRs above the 
Y>* thresholds, Fig. 7 shows that K need not be reduced by much more than by a factor of 10 to 
obtain agreement between theory and all the data points. Further reductions are unnecessary since 
they shift the predicted p*(>ip*) curves below the tp* thresholds. Nevertheless K must be reduced 
by at least a factor of 10. This is illustrated by the dashed red curve in Fig. 7, which shows that 
/5*(>V'*) predicted for K = 0.1xKx enn is incompatible with all the upper limits other than those 
corresponding to #kern =0.25 " and 0.30 " when N max is increased to 2xl0 22 cm -2 . 

The SFR efficiencies can also be decreased by lowering the slope (3 of the Kennicutt-Schmidt 
law, while keeping K=Kxenn- Although Kennicutt (1998a) found /?=1.40±0.15 when he compared 
projected SFRs with the total gas column densities for a sample of 60 galaxies, Gao h Solomon 
(2004) found (3=1 by comparing the SFRs with the H2 column densities of 32 galaxies. Because 
our knowledge of physical conditions in high-z gas is limited, we cannot rule out the possibility 
that (3 = 1 also holds for the atomic gas content of DLAs. Fig. 8 compares the p*(>ip*) curves 
predicted for (3 = 1 and (3 = 0.6 with the standard case, (3=1 .4. With the possible exception of 
DLAs with 6>dia=0.25 ", the curves predicted for (3=1 are inconsistent with the data. Fig. 8 also 
shows consistency with all the data points is achieved if (3 < 0.6. While we cannot rule out such low 
values entirely, the cases (3=1.0 or 1.4 are at least physically motivated (see Elmegreen 2002 and 
Kravtsov 2003). Because we are unaware of any evidence or physical argument to justify values of 
(3 less than 0.6, we focus instead on other mechanisms for reducing SFR efficiencies in DLAs. 

The SFRs would be reduced if the critical surface density for the Kennicutt-Schmidt law 
(see Eq. 1) were higher in DLAs. The physical basis for this follows from the Toomre instability 
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Fig. 8. — Cumulative comoving SFR density versus ip*. Data points as in Fig. 7. Theoretical curves 
correspond to K = Kxenn for slopes (3= 1.4, 1.0, and 0.6. 
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criterion, which is characterized by the parameter Q= Ko g j^G[imu^ j_ (Toomre 1964), where a g 
is the gas velocity dispersion, and k is the epicyclic frequency. The critical H I column density, 
N c ™ t =K<j g I 'ttG /J-rriH , is set by the condition Q=l: gas with N± > Nf lt is unstable to collapse 
and possible subsequent star formation, while gas with N±<Nf lt is stable. In a disk with a flat 
rotation curve n{r) = \^2v rot /r, where r is radial distance from the center of the disk and v rot is 
the rotation speed. Following Kauffmann (1996) and Mo et al. (1998) we assume the DLA disk 
radius r^ a ~0.1r2oo where r2oo is essentially the virial radius of the dark-matter halo enclosing the 
disk and r^ia is denned by the condition N(r ( ii a ) = N m i n =2x 10 20 cm -2 . Adopting the spherical 
collapse model for disk formation, one finds ty o t/ r 200 = WH(z) (Mo et al. 1998) where H{z) is the 
Hubble parameter. The presence of the H(z) term implies that N±(r) at z =3 must be at least 
a factor of 5 higher than at z=0 for gas to form stars. Specifically, N±(r) must exceed Nf lt (r)= 
2.2xl0 21 (r^ a /r) cm" 2 for gas at radius r to be Toomre unstable when z = 3. But the double 
power-law form of f(N,X) rather indicates the gas is stable. For randomly oriented disks, the 
break of the double power-law occurs at Nd=No where No is the maximum value of N±(r): in 
§ 2.1 we found that A^=3.16x 10 21 cm" 2 . Combining Eqs. 5 and 8 with the solution g(N±,X) 

=k3(N±/No)~ 2 , we find that N±(r) = N 1 + 20 [(r 2 — r 2 )/r^ a ] for r < r < r^ia where r 
satisfies the condition N±(ro) = No. Comparison with N^ lt (r) shows that N±(r) < N^ lt (r) for all 
solutions in which N±(r,ua) < N m i n . Therefore, in situ star formation in DLAs may be suppressed 
because they are subcritical disks. 

On the other hand the uncertainties in our model assumptions lead to uncertainties in the 
values of Nf lt . Moreover, uncertainties in the slope a^=dlogf{N,X)/dlogN at N > Nd indicate 
that the disk hypothesis, which predicts «4 = —3, is not well established (see §2.1). Therefore, 
gravitationally unstable gas with > 3.16xl0 21 cm" 2 may be present in DLAs. As a result, while 
cosmological evolution may be an important factor in rendering DLA gas subcritical, additional 
factors may be needed to suppress star formation at column densities up to 10 22 cm" 2 . One such 
factor may be the low molecular content of DLAs: molecules are crucial for star formation since 
emission from rotational transitions of CO and other molecules enables clouds to cool below the 
atomic lower limit of ~ 100 K and to ultimately collapse and form stars. Whereas the median H2 
content in the diffuse neutral gas of the Galaxy, /# 2 ~10 -1 (Tumlinson et al. 2002), Ledoux et al. 
(2003) find that fjj 2 < 10" 6 in DLAs. Interestingly, the median values for fn 2 in the LMC and 
SMC are less than 10 -4 . Tumlinson et al. (2002) argue that the low molecular content of the LMC 
and SMC is likely due to low dust content and relatively strong FUV radiation fields. Since the 
low dust level suppresses H2 formation on grains while the high FUV radiation level increases the 
H2 photodissociation rate, higher values of are required for the gas to be self-shielding against 
photodissociating radiation. Because DLAs have lower dust-to-gas ratios than either the LMC or 
SMC but have similar FUV radiation fields (WPG03), the same mechanisms may by responsible 
for the low molecular content of DLAs. Calculations of photodissociation equilibrium suggest that 
fn 2 would be suppressed at < 10 22 cm" 2 in DLA gas (Tumlinson et al. 2002; Hirashita & 
Ferrara 2005). Krumholz & McKee (2006) and Blitz & Rosolowsky (2006) argue that ip* oc f H2 , 
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which would help to explain the factor of ten or more reduction in the efficiency of star formation 
in DLAs. Consequently, while gas in DLAs with N > N ant < fO 22 cm" 2 is Toomre unstable, it 
may be unable to collapse to form stars. This is in contrast to the Galaxy and most metal-rich 
spirals studied by Martin & Kennicutt (20001) in which the condition N > N CTlt normally results 
in star formation. The explanation (R. C. Kennicutt 2006, priv. comm.) may be the result of a 
coincidence between N CTlt and the column density at which molecular gas becomes self shielding. 
Because of its relatively high dust-to-gas ratio, this is lower in the Galaxy (and most metal-rich 
spirals) than in metal-poor objects such as the LMC, SMC (Tumlinson et al. 2002), and in DLAs. 
This may help to explain why the Toomre criterion alone appears to be a sufficient condition for 
the onset of star formation in most metal-rich spirals (Martin & Kennicutt 2001), but not in DLAs. 

A further mechanism for cloud stabilization was suggested by Ferguson et al. (1998) who 
argued that flaring of the H I disk in the outer regions of spiral galaxies leads to suppression of star 
formation: for a given gas surface density the increased scale-height reduces the volume density, 
which helps to stabilize gas against Jeans instability. Schaye (2004) argued that disks would be 
Toomre stable provided the gas is a warm (T ~ 10 4 K) neutral medium and unstable when the gas 
is a cold (T ~ 100K) neutral medium. But independent lines of evidence (e.g. WGP03; Howk et 
al. 2005) suggest that half of the known sample of DLAs is comprised of cold (T ~ 100K) gas with 
column densities exceeding the (5— 20)xl0 20 cm" 2 instability threshold. While this mechanism 
needs to be investigated further, current evidence suggests that phase transitions alone are not 
sufficient to trigger star formation in DLAs. There may be other mechanisms that do suppress star 
formation, but further discussion is beyond the scope of this paper. 



5.2. Implications for Metal Production and Heat Input: the LBG Connection 

The UDF results place a conservative upper limit on p\ contributed by DLAs with N > N m i n , 
which is given by p* < 10~ 2 ' 4 M Q yr _1 Mpc -3 for z ~ 3 (see curve with K=0.1xKK enn in Fig. 7). 
The new limit on p\ has several consequences. First, since the metal production rate is proportional 
to p\, one can compute the DLA metalicity at a given redshift from the previous history of p*(z). 
Nagamine et al. (2004b) used their numerical simulations, which incorporated the Kennicutt- 
Schmidt law, to predict p\ = 10" 1 M Q yr" 1 Mpc -3 and the metallicity [M/H]=— 0.5 at z=3. If we 
scale their result to our upper limit on p\, their chemical evolution model predicts [M/H] < —1.9. 
By comparison the average DLA metallicity at z ~ 3 is given by [M/H]=— 1.40±0.07 (Prochaska et 
al. 2003). Therefore, our limit on p\ has shifted the problem of metal overproduction by a factor 
of ten in DLAs (see WGP03) to one of underproduction a factor of three. 

Second, the new limit on p\ implies a gas heating rate significantly lower than indicated by 
the measured cooling rate. The [C II] 158 pm cooling rate per unit comoving volume inferred from 
C II* absorption in DLAs (Wolfe et al. 2004) is given by 



C = Q g Pcrit< 4 >/fJ>m H 



(9) 
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where fl g is the mass per unit comoving volume of neutral gas in DLAs, < t c > is the average [C 
II] 158 /im cooling rate per H atom, and //m# is the average mass per particle. Assuming p=1.3, 
our data indicate C= (2±0.5)xl0 38 ergs s _1 Mpc~ 3 in the redshift interval z= [2.5,3.5]. Here we 
have ignored half the DLA population with £ c < 10 -271 ergs s _1 H _1 because it is likely heated 
by background radiation alone (Wolfe 2005). 

To balance cooling, Wolfe et al. (2003a [hereafter WPG03]) suggested the grain photoelectric 
effect as a plausible heating mechanism. The comoving heating rate is given by 



where §{L v )dL v is the comoving density of DLAs with FUV luminosities in the interval (L u , L v +dL u ) 
The heating rate per DLA is given by 



where T(r) is the heating rate per H atom, n(r) is the density of H atoms at displacement vector 
r, and the integral extends over the DLA volume. From the models of Bakes & Tielens (1994) and 
Weingartner & Draine (2001), WPG03 find that r=10" 5 KeJ^(r) where k is the dust-to-gas ratio, e 
is the heating efficiency, and J u (r) is the FUV mean intensity. Assuming DLAs are uniform disks 
with radius R and scale-heights h, and adopting the solution for J u (WPG03) one finds 



where < N > is the average H I column density of the DLA sample. Because p*=1.2f>x 
10~ 28 J" 4>{L u )L v dL u (Madau et al. 1998), the comoving heating rate TL is proportional to p'*. 
Assuming parameters adopted by Wolfe et al. (2004), we find that our upper limit on p* implies 
TL < 4xl0 37 ergs s _1 Mpc~ 3 . Therefore, in situ star formation in DLAs is unlikely to account for 
the energy input needed to balance cooling. 

Consequently the new limits on p'* require an external source of metals and heat input for 
DLAs with i c > 10 -27 ' 1 ergs s _1 H _1 . This idea was first discussed by WGP03 who suggested that 
centrally located galaxy bulges could heat the absorbing gas. Wolfe (2005) then suggested LBGs 
as the most natural bulge sources. In this scenario the FUV radiation originates in the LBGs, 
is attenuated by their high dust content, and then propagates unattenuated throughout DLAs 
because of their low dust content. The LBGs heat the surrounding DLA gas, with the same FUV 
radiation observed by us to have p'*= 10~ L7 M Q yr -1 Mpc~ 3 at z=3 (e.g. Giavalisco et al.2004). 
From Eq. 12, which also applies for central point sources (WGP03), we find the resulting LBG 
heating rate to be W=(3±2)xl0 38 erg s _1 Mpc -3 , where the large systematic errors stem from 
uncertainties in the dust composition, photoelectric heating efficiency, geometry, etc. Nevertheless 




(10) 




(11) 




(12) 



-26- 



this estimate suggests that embedded LBGs can explain the heating rate inferred for DLAs with 
4 > 1CT 27 - 1 ergs s" 1 H" 1 . 

The physical association of high- z DLAs with LBGs is supported by the following phenomena. 
First, Cooke et al. (2006) recently measured the 3-D cross correlation function between DLAs and 
LBGs and found good agreement between the DLA-LBG cross-correlation amplitude and the LBG- 
LBG autocorrelation amplitude. This implies a similarity in the dark-matter masses of the two 
populations; i.e., a strong overlap between DLAs and LBGs. Further evidence for this association 
comes from the identification of at least one high-z DLA with an LBG (M0ller et al. 2002). In 
this object the FUV surface-brightness of the compact emitting regions is two orders of magnitude 
higher than predicted by the Kennicutt-Schmidt law applied to the H I column-density detected 
in absorption. This suggests that star formation originates in compact, presumably dusty, regions 
rarely detected in absorption, which are surrounded by the lower column-density gas that gives 
rise to the DLA. This configuration brings to mind the "tadpole galaxies" discussed in § 3.2 and 
indicates that their star forming regions may also be embedded in DLA gas. The low area covering 
factor and high column density implies that star formation in DLAs occurs in dense molecular 
regions (see Zwaan & Prochaska 2006), which are replenished with gas from the surrounding DLA. 
By contrast, the DLA gas may be chemically enriched by P-cygni winds emanating from the star- 
forming LBGs (Steidel et al. 2003). Other sources of metals may be pre-enriched gas or stars 
accreted by the galaxy hosting the DLA. The latter processes may be crucial for DLAs with i c < 
10 -27 - 1 erg s _1 H _1 , which may not contain centrally located LBGs. 



6. SUMMARY and CONCLUDING REMARKS 

We searched for in situ star formation in the neutral gas comprising DLAs by looking for 
low surface-brightness emission from spatially extended objects in the Hubble Ultra Deep Field 
F606W image. The search was designed to detect objects in the redshift interval z = [2.5, 3.5] with 
linear diameters ranging between 1.9 kpc and 31 kpc, which encompasses most model predictions 
published so far. The search was optimized using a smoothing kernel to match the predicted 
size of the DLAs, ranging between an angular diameter of #dia= 0.25" and 4.0". The results are 
summarized as follows: 

• After eliminating compact objects with high surface brightness, we found the number of 
extended sources to decrease rapidly with increasing #kern: the implication is that extended sources 
of low surface-brightness emission rarely occur in the UDF (see Fig. 3). Because photometric 
redshifts are unavailable, these findings lead to upper limits on the comoving density of objects 
within the UDF with z= [2.5,3.5]. 

• The upper limits on comoving density lead to upper limits on comoving SFR density, which 
are between factors of 30 to 100 lower than predicted by the Kennicutt-Schmidt law. As a result, 
the projected SFRs per unit area must be at least a factor 10 lower than the rates predicted by the 
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distribution of neutral-gas column densities in DLAs. The upper limits are quite general as they 
apply to objects with angular diameters encompassing the sizes of DLAs in ACDM and all other 
models suggested so far. 

• We suggest that lower SFR efficiencies may partly result from an increase with redshift of 
the critical surface density for Toomre instability. Such an increase is due to cosmological evolution 
alone. This effect causes the DLA gas to be Toomre stable, which could remove the discrepancy 
between the comoving SFR density predicted by the Kennicutt-Schmidt law and our empirical 
upper limits. However, the uncertainties are such that gravitational instabilities may be present in 
DLAs. In that case we suggest that the low molecular content of the gravitationally bound clouds 
with N < 10 22 cm -2 prevents further cooling and thus inhibits star formation. The low molecular 
content of the gas may be related to the minimum column density required for molecule formation 
(Krumholz & McKee 2005; Blitz & Rosolowsky 2006) since molecule formation in objects with 
low dust content and high FUV radiation fields such as found in DLAs could require values of 
./V exceeding 10 22 cm' 2 for shielding against photo-dissociating radiation. The low area covering 
factor of such high column densities may be the reason for the low molecular content measured in 
DLAs. 

• The upper limits on comoving SFR density reduce the rate of metal production predicted in 
DLAs. The resultant upper limits on metallicity at z ~ 3 are lower than the observed metallicities 
by a factor of three. While this removes the previous problem of metal overproduction (by a factor 
of ten) in DLAs, it leads to a new problem of metal underproduction. 

• The upper limits on comoving SFR density result in a comoving grain photoelectric heating 
rate that is significantly lower than the comoving cooling rate inferred from [C II] 158 /xm emission 
rates for about half the DLA sample. Therefore, external sources of heating are required. We 
argue that the observed comoving SFRs of LBGs with R < 27 and z= [2.5,3.5] indicate a heat input 
sufficient to balance cooling. 

• We suggest a scenario in which half the DLA population with significant [C II] 158 /um cooling 
is powered by centrally located, compact LBGs. We speculate that DLA metal enrichment may be 
due to P-cygni outflow from the metal-enriched LBGs. The half of the DLA population without 
significant [C II] 158 /um cooling is presumably heated by background radiation alone. Without an 
LBG, the source of metals for this DLA population is unclear. 

Two possible caveats come to mind. First, since our limits do not apply to emission from DLAs 
with #dia < 0.25 it is possible that for most DLAs at z ~ 3, d^ia is less than 1.9 kpc and the 
SFRs follow the Kennicutt-Schmidt relation. However, as discussed in § 3.1 all models suggested 
so far predict that the bulk of DLAs will have diameters larger than 1.9 kpc. 

The second caveat stems from the possibility that in situ star formation at the Kennicutt- 
Schmidt rate occurs in DLA gas associated with high surface-brightness objects, since such objects 
were automatically excluded from the survey. We tested this hypothesis by searching for low 
surface-brightness emission surrounding high surface-brightness objects. To enhance the sensitivity 
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of our search we smoothed the regions in the F606W image surrounding selected targets using 
kernels with 0.5 " < #kem < 2.0 ". In all cases the convolved images failed to show evidence for 
outlying emission with fluxes comparable to the central object. But it is difficult to draw conclusions 
from these data since the redshifts of all but two galaxies, at z = 2.593 (Szokoly et al. 2004) and 
z = 3.797 (Vanzella et al. 2006), respectively, are unknown. We found these two galaxies to consist 
of a compact core (FWHM = 0.3 ") surrounded by diffusion emission, which exhibited maximum 
contrast with respect to sky when #kern = 0.5 ". In the case of the z=2.593 galaxy the magnitude 
of the diffuse emission, V=30, while the central core has V = 25; i.e., the regions surrounding this 
star-forming galaxy contributes less than 1 % of the total SFR. Because p\ due to compact star 
forming galaxies is comparable to p\ predicted by the Kennicutt- Schmidt relation for DLAs, we 
would conclude that neutral gas surrounding star- forming objects does not form stars according to 
the Kennicutt-Schmidt relation if this galaxy is typical of its class. 

Finally we compare our empirical limits with previous analyses that apply the local Kennicutt- 
Schmidt law for evaluating p\ in DLAs. First, Lanzetta et al. (2002) measured the frequency 
distribution of projected SFR per unit area, h(ip*), for galaxies in the HDF in the redshift interval 
2;= [0,10]. Their measurements at z ~ 3 were sensitive to > 10 -1 ' 5 M Q yr _1 kpc~ 2 . To determine 
h(ijj*) at lower values of ip*, they applied the Kennicutt-Schmidt law to the H I column-density 
distribution function f(N,X) by using the relation h(ip*)dip*=(Ho/c)f(N,X)dN. They then used 
Eq. 6 in the spherically symmetric limit to determine p\. These authors concluded that at z k, 3 
the bulk of SFR density resides in regions of 0.1 M yr _1 kpc~ 2 , where the majority of LBG's 
are found. They argued that DLAs of N < 5 x 10 21 cm~ 2 contribute an insignificant fraction 
of the global SFR density. This is inconsistent with the current form of f(N, X) which when 
extrapolated to iV » 2xl0 22 cm~ 2 , to include HSB objects such as LBGs, predicts that more 
than 50 % of p\ is contributed by DLAs with N < 5 x 10 21 cm~ 2 , even when i^=K^ enn /2.8 as 
suggested by Lanzetta et al. (2002). The discrepancy stems from the different expressions used to 
compute f(N, X): Lanzetta et al. (2002) used the then available expression for f(N, X) with slope, 
03=— 1.4, that is flatter than 03=— 2.0 derived for the more accurate SDSS function (PHW05) used 
here. Because we have shown that the DLA contribution to p\ cannot be this large, the newer data 
do not support the suggestion of Lanzetta et al. (2002) that SFRs predicted by application of the 
Kennicutt-Schmidt law to DLAs overlap those derived directly from optical emission. On the other 
hand, their analysis is consistent with our interpretation that LBGs provide the bulk of the star 
formation at z ~ 3. Second, Hopkins et al. (2005) combined an updated version of f(N, X) with 
the Kennicutt-Schmidt relation to predict p\ in the spherically symmetric limit for the redshift 
interval z= [2.5,3.5]. These authors concluded that DLAs contribute roughly 50% of the observed 
SFR density at z ~ 3. Again, the predicted p\ is significantly higher than our empirical upper 
limits, and thus their models are ruled out by the UDF data. The new upper limits on p\ are also 
inconsistent with the predictions of the uniform-disk model of WPG03 and WGP03. Their models 
do not make use of the Kennicutt- Schmidt law, but rather use the strength of C II* absorption to 
infer the level of in situ star formation throughout DLAs modeled as uniform disks. Because the 
predicted average projected SFR, <^*>=10 -2 M Q yr _1 kpc -2 , is above our detection threshold, 
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the factor of 100 disagreement between the predicted p* and our upper limits rules out this model. 
On the other hand, the new upper limits are consistent with the "bulge model" suggested by 
WGP03. In this model, star formation does not occur in the gas detected in absorption, but rather 
arises in a centrally located compact region, which heats the surrounding disk gas by emitting 
FUV radiation. The rough agreement between the comoving heating rate supplied by LBGs and 
the comoving cooling rate of DLAs suggests the bulges are in fact LBGs. 

We wish to thank Eric Gawiser, Kim Griest, Andrei Kravtsov, and Jason X. Prochaska for 
valuable discussions. AMW was partially supported by NSF grant AST 03-07824 and H-WC was 
partially supported by NASA grant NNG06GC36G. 
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